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慧眼 -HXMT 观 测 到 Don 双星 系统 中 的 W-K 
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摘要 ”通过 几 十 年 的 观测 


究 , 黑洞 X 射 线 双 星 (X-Ray Binary, XRB) 部 分 特征 被 揭示 . 然而 ， 


不 确定 . 黑洞 XRB 功 率 密度 谱 的 截 
系 ) 可 以 限制 吸 积 盘 结构 . 利 
据 , 对 黑洞 XRB 的 W-K 关 系 进行 了 厂 
在 MAXI J1535-571 之 中 存在 截 
功率 密度 谱 来 自 质 
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量 吸 积 率 的 扰动 传播 , 可 以 推 


期 振荡 (Quasi Periodic Oscillation, QPO) 的 相关 
民 -HXMT (Hard X-ray Modulation Telescope) 观 测 到 的 5 个 黑洞 XRB 的 数 
在 慧眼 -HXMT 观 测 的 3 个 探测 器 能 段 中 W-K 关 系 成 立 . 此 外 
断 频 率 和 吸 积 盘 内 半径 的 相关 关系 , 这 和 截断 
则 吸 积 盘 内 半径 接近 最 内 


科学 院 高 能 物理 研究 所 北京 100049) 


及 积 盘 结构 尚 
性 质 (W-K 关 


的 吸 积 盘 结 构 一 致 . 如 果 观 测 到 的 
区 稳定 轨道 , 此 黑洞 可 能 是 高 自 旋 
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关键 词 XH: 双星 , 相对 论 过 程 , 吸 积 盘 , 辐射 机 制 : 普通 
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1 引言 

X 射 线 双 星系 统 (X-Ray Binary, XRB) 通 常 是 
由 一 个 致密 天 体 和 一 个 伴星 组 成 , 根据 伴星 质 
量 的 大 小 可 以 分 为 小 质量 X 射 线 双 星 (Low-Mass 


X-ray Binary, LMXB) 和 大 质量 X 射 线 双 星 (High- 
Mass X-ray Binary, HMXB) 系 统 . 根据 XRB 中 的 
致密 天 体 可 以 分 为 黑洞 Xx 射线 双星 (Black Hole X- 


度 Lx < 10°! erg- s71), 有 时 候 由 于 吸 积 流 的 热 不 
稳定 性 导致 源 产生 亮 的 爆发 现象 . 

黑洞 暂 现 源 爆 发 的 演化 随 着 光度 的 变化 而 表 
现 出 不 同 的 谱 态 (spectral state)?), 其 能 谱 成 分 很 
复杂 , 但 辐射 能 谱 通 常 简单 地 划分 为 两 部 分 , 来 
自 吸 积 盘 的 热 谱 成 分 和 热 内 流 ( 冕 ) 的 非 热 谱 成 分 . 
在 爆发 的 开始 阶段 , 源 通常 处 于 非 热 辐射 成 分 主 


Ray Binary, BHB)、 中 子 星 X 射 线 双 星 等 . 由 于 
伴星 物质 被 中 心 致密 天 体 吸 积 , 其 引力 势能 通 
过 物质 的 相互 作用 转化 为 辐射 能 , 主要 辐射 能 
在 X 射 线 波段 . 洛 希 激 的 吸 积 会 在 中 心 致密 天 体 
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时 的 低 硬 态 (Low Hard state, LHS), fet 7% 48 
By ~ 1-1.7, 能 谱 在 100-200 keV 处 截断 . 随 着 物质 
吸 积 率 的 增加 , 热 谱 成 分 逐渐 增强 , 源 演化 进入 中 
间 态 (Intermediate state, IMS). 基于 能 谱 和 时 变性 
质 , 我 们 可 以 进一步 将 IMS 划 分 为 硬 中 间 态 (Hard 
Intermediate state, HIMS) 和 软 中 间 态 (Soft Inter- 
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mediate state，SIMS)B. 随 着 光度 的 进一步 增 
WM, 当 热 成 分 主导 源 的 辐射 时 , 源 进入 了 高 软 态 
(High Soft state, HSS), 在 高 软 态 源 的 7 一 般 大 于 2， 
一 般 认 为 吸 积 盘 内 半径 处 于 最 内 稳定 圆 形 轨 
道 (Innermost Stable Circular Orbit, ISCO). 在 


学 R 
型 , 现在 有 越 来 越 多 的 观测 证 据 支持 C 型 QPO 来 自 


热 内 流 相 对 论 进 动 模型 中 的 L-T 进 动 (Lense- 
Thirring procession)l13; 15716), L-THH Xt W HE By He 
型 最 初 被 Stella 等 ("用 来 解释 中 子 星 LMXB 中 
的 QPO 现 象 , 在 此 之 后 由 Ingram 等 L 引 将 其 应 用 在 


爆发 结束 时 , 源 经 历 中 间 态 返回 低 硬 态 , 或 者 宁 
静态 辐 . 黑洞 LMXB 的 一 次 完整 的 爆发 过 程 , 也 
可 以 用 硬度 -强度 图 (Hardness-Intensity Diagram， 
HID) 进 行 体现 , 其 中 硬度 是 高 能 段 光 子 计数 率 和 
低能 段 光 子 计数 率 的 比率 . 对 于 典型 黑洞 LMXB 的 
一 次 爆发 过 程 , 源 的 演化 在 HID 图 上 面 呈现 出 一 个 
逆 时 针 的 q 形 状 . 基于 爆发 在 HID 的 位 置 , 也 可 以 将 
爆发 按照 前 述 分 成 4 个 谱 态 : 低 硬 态 、 硬 中 间 态 、 
软 中 间 态 和 高 软 态 . 同时 也 展示 了 谱 态 在 爆发 过 程 
中 的 演化 , 理论 研究 倾向 于 认为 这 种 变化 是 由 于 
吸 积 盘 的 截断 半径 向 内 移动 和 吸 积 率 的 增加 引起 
gl, 


黑洞 X 射 线 双 星 表 现 出 广泛 的 长 时 标的 能 谱 
状态 变化 和 短 时 标 时 变现 象 . 对 黑洞 暂 现 源 的 时 
变 研究 发 现 , 准 周期 振荡 (Quasi Periodic Oscilla- 


黑洞 XRB 之 中 05719. 在 此 模型 中 , 处 于 吸 积 盘 的 截 
断 半 径 之 内 的 热 内 流 受 黑洞 强 引 力 场 和 自 旋 的 影 
响 产 生 相 对 论 进 动 , 即 观测 到 的 低频 QPO 现 象 . 
通过 对 时 变性 质 的 研究 发 现 , 黑洞 和 中 子 星 
LMXB 的 QPO 中 心 频率 与 噪声 特征 截断 频率 
在 对 数 坐 标 图 上 呈 线 性 正 相 关 ， 且 两 类 天 体 
遵循 相同 的 规律 这 一 关系 被 Wijnands 和 van 
der Klis!) HN, 简称 为 W-K 关 系 . W-K 关 系 表 
明 , 在 所 有 这 些 不 同类 型 的 源 中 , 带 限 噪 声 和 
QPO 是 由 相同 的 物理 机 制 产生 的 . 它 很 可 能 
与 内 部 的 吸 积 流 的 不 稳定 性 有 关 , 而 与 致密 天 体 的 
质量 与 磁场 无 关 . 这 表明 黑洞 和 中 子 星 附近 吸 积 流 
可 能 遵循 相似 的 物理 规律 , 同时 也 为 研究 致密 天 体 
吸 积 流 的 运动 规律 提供 了 途径 . 
a HR-HXMT (Hard X-ray Modulation Tele- 


ms 


i 普遍 存在 于 黑洞 LMXB， 功率 密 
tion, QPO) 普 裔 存 在 黑洞 LMXB. QPO 在 ea scope) 卫 星 自 2017 年 发 射 以 来 , 已 成 功 对 黑洞 暂 
度 谱 (Power Density Spectrum, PDS) ERI NÆ 和 Wats i D 
i Saye a eee AS 现 源 进行 了 广泛 观测 . 本 文 使 用 了 慧眼 -HXMT 在 
的 峰值 , 其 他 部 分 被 称 为 宽带 噪声 . QPO 分 为 低频 eA ee we mee Tia 
29 ee 2017 一 2020 年 间 观 测 的 5 个 典型 黑洞 暂 现 源 的 数据 ， 
QPO 和 高 频 QPO, 低频 QPO 又 可 以 通过 PDS 图 上 Ea 
Fn ; > 研究 了 其 时 变性 质 . 慧眼 -HXMT 的 宽 能 带 特征 (1- 
状 、 均 方 根 (Root Mean Square, RMS)、 品 质 因 CARRS AA Set a RS 
A Pope A RSH anys 250 keV) 可 以 让 我 们 在 更 宽 的 能 量 范围 研究 黑洞 
子 和 时 间 延 迟 等 性 质 分 为 A、B 和 C 这 3 类 [6-9), 我 们 oe ie sorta 
= Se bee tea j F 暂 现 源 的 能 谱 和 时 变性 质 . 我 们 的 研究 将 W-K 关 系 
般 用 rms 直 接 代表 “相对 RMS”. 其 中 C 型 QPO 伴 拓展 到 120 keV 的 能 段 并 对 截断 盘 的 模型 做 出 限 
随 着 平台 形状 的 低频 噪声 成 分 ( 称 之 为 带 限 品 E eee acer 


FA), 这 一 类 型 QPO 出 现 于 硬 态 和 硬 中 间 态 . 通 
过 比较 黑洞 和 中 子 星 LMXB 的 QPO 现 象 人 们 认 
Wy, 它们 可 能 具有 相似 的 物理 起 源 [6 10121, 对 于 
本 文 没 有 详细 涉及 的 其 他 类 型 QPO, 可 以 参考 
Motta 等 中 、Wijnands 等 6、Psaltis 等 0、Belloni 
等 上 的 工作 . 

对 于 X 射 线 双 星 中 的 低频 准 周 期 振荡 (Low 
Frequency Quasi Periodic Oscillation, LFQPO) 进 
行 了 广泛 研究 后 , 提出 了 多 个 模型 来 解释 它们 的 起 
RUS 14) | 这 些 模型 一 般 分 为 吸 积 的 几何 和 物理 模 


制 |， 
2 数据 处 理 


慧眼 -HXMT 是 中 国 在 2017 年 6 月 15 日 成 功 发 
射 的 第 一 颗 X 射 线 天 文 望远镜 Ha. 慧眼 -HXMT 由 3 
个 准 直 望远镜 组 成 : 高 能 X 射 线 望 远 镜 (HE)D9)、 
中 能 XxX 射线 望远镜 (ME)R9q 和 低能 Xx 射线 望 远 
镜 (LE)P， 具 有 能 段 宽 、 有 效 面积 大 、 时 间 分 辩 
率 高 等 特点 , 适用 于 黑洞 时 变 和 能 谱 性 质 研究 . 臣 
眼 -HXMT 基 本 参数 可 见 表 1. 
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马 宾 远 等 : 慧眼 -HXMT 观 测 到 的 


星系 统 中 的 W-K 关 系 


黑洞 双 


表 1 慧眼 -HXMT 3 个 探测 器 的 参数 
Table 1 Parameters of the three payloads of Insight-HXMT 


Detector LE ME HE 
Geometrical area/cm? 384 952 5100 
Energy band/keV 1-12 8-35 20-350 
Time resolution/s 1000 240 4 
Energy resolution 2.5% @6keV 15% @ 20keV < 17% @ 60 keV 


为 了 探究 黑洞 双星 的 时 变性 质 , 我 们 使 用 了 
慧眼 -HXMT 对 5 个 黑洞 双星 (MAXI J1535-571、GX 
339-4, MAXI J1348-630, MAXI J1631-479 和 GRS 
1915+105) 的 观测 数据 , 根据 已 有 的 工作 B2726] 对 
源 的 数据 分 析 和 谱 态 划分 , 我 们 仅 选 取 了 有 明 
确 C 型 QPO 的 LHS 和 SIMS 的 观测 号 ( 表 2-4), 表 中 
从 左 到 右 分 别 是 : 观测 号 (ObsID), 时 间 (MJD)， 
低频 截断 频率 (ww,),， 高 频 截断 频率 (ww)，QPO 频 
率 (zopo)， 其 中 有 些 观 测 号 在 某 一 能 段 上 有 和 缺 
失 . 我 们 用 慧眼 -HXMT 数 据 分 析 软 件 (HXMTDAS 
v2.041) 对 表 2-4 中 所 有 数据 进行 分 析 , 各 项 参数 使 
软件 默认 值 : (1) 指 向 偏 移 角 (offset angle of the 
pointing direction) < 0.04°; (2) 仰 角 (earth eleva- 
tion of FOV center position, ELV) > 10°; (3) 地 
磁 截 止 刚 度 (geomagnetic cut-off rigidity, COR) > 
8 GV. 

对 高 能 天 体 辐射 的 时 变 分 析 , 通常 采用 传 里 叶 
分 析 方 法 计算 光 变 曲线 的 功率 密度 、 时 延 和 相关 
系数 等 频 域 特征 量 0 271, 在 这 里 我 们 使 用 powspec 
软件 进行 功率 密度 谱 的 计算 , 设 定 PDS 由 128 s 的 
时 间 间 隔 (interval) 产 生 ， 时 间 并 (bin) 为 1/128 s, 
即 Nyquist 频 率 为 64 Hz. 然后 对 全 部 间隔 的 PDS 进 
行 平均 , 得 到 每 个 观测 号 的 PDS, FF EVRA TAPS 
声 后 进行 Miyamoto 归 一 化 Pa. 死 时 间 效 应 在 LE、 
HE 可 以 忽略 . ME 的 死 时 间 处 理 方 法 参考 Huang 等 
ARTERA 为 增加 PDS 的 信 噪 比 , 观测 时 间 小 
于 640 s 的 数据 被 忽略 , 同时 , QPO 频 率 相近 (中 心 


*http://hxmtweb.ihep.ac.cn/software.jhtml 


频率 差异 小 于 QPO 的 半 高 全 宽 ) 且 时 间 相 近 的 观测 
被 合并 计算 PDS. 在 我 们 所 选择 HXMT 观 测 源 中 ， 
MAXI J1535-571、GX 339-4, MAXI J1348-630、 
MAXI J1631-479 和 GRS 1915+10547 GA 10, 2, 8, 
3, 3 个 有 效 观 测 号 ( 表 2-4). 

如 引言 部 分 所 述 , LHS 和 HIMS 的 PDS 由 宽带 
噪声 和 QPO 构 成 . 我 们 选择 使 用 XSPEC v12.12.12, 
按照 Wijnands 等 回 拟 合 PDS 的 方式 使 用 Lorentz 
模型 拟 合 QPO 及 其 谐 频 成 分 、 使 用 bknpow 或 
者 bkn2pow 两 种 模型 拟 合 噪声 成 分 , 得 到 QPO 频 
率 和 宽带 噪声 截断 频率 . 拟 合 频率 范围 是 
0.01-64 Hz. 图 1 展示 了 两 个 典型 例子 , 上 图 是 观测 
号 P0114535005 的 ME 观测 的 PDS, 在 其 宽带 噪声 上 
可 以 看 到 一 个 明显 的 截断 , 我 们 用 bknpow 进 行 拟 
合 , 得 到 截断 频率 vw; 下 图 是 观测 号 P0114535005 
的 ME 观测 的 PDS, 用 bkn2pow 进 行 拟 合 , 可 以 得 到 
两 个 截断 频率 和 wh. XSPEC 的 “error” 命 令 给 出 
各 个 参数 的 误差 , 默认 的 置信 和 度 为 90%, 如 表 2-4 中 
截断 频率 和 QPO 频 率 参 数 的 误差 . 对 于 本 文 处 理 
的 观测 号 , 低 硬 态 的 ww 一 般 为 0.04-0.4 Hz, 而 在 
人 硬 中 间 态 可 以 达到 和 ~1-2 Hz. Lorentz 模 型 的 中 心 
频率 就 是 QPO 的 频率 vopo. XSPEC 在 拟 合 完 成 
后 还 可 以 给 出 拟 合 的 约 化 x? 值 (反映 了 拟 合 的 优 
劣 , 越 接近 1 则 拟 合 越 好 ), 我 们 的 拟 合 中 约 化 x? 值 
小 于 1.5 (自由 度 为 186), 大 多 数 在 1.2 左 右 ， 可 见 
拟 合 效果 较 好 . 我 们 使 用 XSPEC 的 命令 “fux” 估 
计 2-120 keV 的 不 同 能 量 段 的 低频 噪声 成 分 的 rms. 
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#2 慧眼 -HXMT 的 LE 观测 到 的 5 个 BHB 在 2-10 keV 能 段 的 PDS 拟 合 得 到 的 结果 
Table 2 The PDS fitting results of five BHBs observed by Insight-HXMT/LE in 2-10 keV 


ObsID MJD Vp,LE/ Hz ZhLE/Hz Vepo,LE/Hz 


MAXI J1535-571 
P011453500108 58002.734 0.26 +0.02 4.344 


| 
iT 
> 
N 
A 


1.31 + 0.02 
P011453500144 58008.444 0.43 0.05 9.02 Æ 0.65 2.61 + 0.01 
P011453500145 58008.583 0.30 0.02 10.31 0.97 2.60 0.01 


P011453500301 58011.200 0.36 +0.05 8.09 0.68 2.12 + 0.01 
P011453500401 58012.260 0.3440.03 8.454 1.28 2.75 + 0.01 
P011453500501 58013.255 0.45 0.03 14.53 3.17 3.35 0.01 


P011453500601 58014.117 0.63 0.04 9.27 + 1.86 3.30 0.01 


P011453500901 58015.974 1.48 0.05 xx? 8.74 0.03 
P011453500902 58017.250 1.73 0.05 六 站 9.00 + 0.03 
P011453500903 58017.389 1.45 0.05 ex 7.74 + 0.05 
GX 339-4 
P0304024035 59295.758 0.12 0.13 ex 0.64 + 0.01 
P0304024036 59297.017 0.11 0.03 1.3340.12 0.93 + 0.01 
MAXI J1348-630 
P0214002002 58511.314 0.02+0.01 5.49 + 0.28 0.39 + 0.01 
P0214002003 58512.308 0.1140.01 5.96 + 0.27 0.54 0.00 


P0214002004: 01-08? 58513.303 0.16 +0.01 6.26 0.19 0.71 0.06 
P0214002004: 09-19 58514.4111 0.12+0.01 7.08 0.25 0.86 + 0.02 


P0214002005 58516.818 0.22+0.01 12.190.511 1.25 0.01 
P0214002006 58518.012 0.36 +0.04 4.57 0.31 2.81 + 0.05 
P0214002007 58519.145 0.56 +0.03 2.78 + 0.09 3.82 + 0.04 
P0214002008 58519.804 0.44+0.16 2.97 + 0.27 3.37 0.12 


MAXI J1631-479 


P0214003002 58525.973 0.84 + 0.06 ex 4.85 + 0.05 

P0214003003 58527.1 1.19+0.40 13.7541.70 6.21 0.09 

P0214003008 58534.787 4.73 + 0.36 ex 6.78 + 0.35 
GRS 1915+105 

P0101330006 58216.956 1.3140.13 9.60 1.21 3.67 0.02 

P0101330007 58218.216 0.52 0.18 ex 4.20 + 0.04 

P0101330009 98220.272 1.15 + 0.23 ex 3.39 + 0.02 


a ** represent the PDS does not have this noise component in fitting. 


> The number behind the colon represents the exposure ID in this observation. 
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表 3 慧 眼 -HXMT 的 ME 观测 到 的 5 个 BHB 在 10-20 keV 能 段 的 PDS 拟 合 得 到 的 结果 
Table 3 The PDS fitting results of five BHBs observed by Insight-HXMT/ME in 10-20 keV 


ObsID MJD Vp, MBE / Hz Ln,ME/Hz VQPO,ME/Hz 


MAXI J1535-571 
P011453500108 58002.734 0.26+0.01 5.70 土 0.30 1.30 土 0.02 
P011453500144 58008.444 0.6140.08 14.24 土 1.96 2.58 土 0.01 
P011453500145 58008.583 0.53 土 0.04 12.07 土 1.36 2.62 土 0.01 
P011453500301 58011.200 ”0.28 土 0.03 6.23 + 0.27 2.07 + 0.00 
P011453500401 58012.260 0.36+0.03 9.02 40.37 2.76 + 0.00 
P011453500501 58013.255 0.4340.05 7.39 0.35 3.37 + 0.01 
P011453500601 58014.117 0.82+0.06 6.86 + 0.26 3.35 + 0.01 


P011453500901 58015.974 1.58 + 0.04 xx? 9.05 0.01 
P011453500902 58017.250 1.70 0.05 ex 9.28 + 0.01 
P011453500903 58017.389 ex 六 六 六 六 
GX 339-4 
P0304024002 59264.083 0.42 + 0.03 ex 1.42 + 0.06 
P0304024036 59297.017 0.05 + 0.04 六 六 0.95 士 0.01 
MAXI J1348-630 
P0214002002 58511.314 0.05+0.01 6.73 0.39 0.39 + 0.01 
P0214002003 58512.308 0.1340.01 7.364 0.28 0.57 + 0.02 
P0214002004: 01-08? 58513.303 0.13+0.01 7.73 0.18 0.74 + 0.01 
P0214002004: 09-19 58514.411 0.17+0.01 六 六 0.87 + 0.01 
P0214002005 58516.818 0.22+0.02 12.19+0.52 1.26 0.01 
P0214002006 58518.012 0.35 土 0.05 3.57 0.27 1.39 + 0.02 
P0214002007 58519.145 0.72 土 0.12 5.73 + 0.20 3.69 0.01 
P0214002008 58519.804 0.150.111 5.72 + 0.26 3.4 + 0.01 


MAXI J1631-479 


P0214003002 58525.973 1.387 + 0.37 ex 4.99 + 0.03 
P0214003003 58527.1 1.80 +0.54 10.24+0.80 6.20+0.05 
P0214003008 58534.787 1.70+ 0.97 ex 6.63 + 0.08 


GRS 1915+105 


P0101330006 58216.956 0.45 + 0.03 ex 3.72 + 0.03 
P0101330007 58218.216 0.65 + 0.04 ex 4.32 + 0.13 
P0101330009 58220.272 0.27 + 0.06 ex 3.45 + 0.03 


a ** represent the PDS does not have this noise component in fitting. 


> The number behind the colon represents the exposure ID in this observation. 
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表 4 慧眼 -HXMT 的 HE 观测 到 的 5 个 BHB 在 20-120 keV 能 段 的 PDS 拟 合 得 到 的 结果 
Table 4 The PDS fitting results of 5 BHBs observed by Insight-HXMT/HE in 20-120 keV 


ObsID MJD Vp HE / Hz Mp HE / Hz Vepo,HE/ Hz 


MAXI J1535-571 
P011453500108 58002.734 0.49 + 0.22 xx? 1.50 0.03 
P011453500144 58008.444 0.51 0.18 7.87 + 0.51 2.56 + 0.01 
P011453500145 58008.583 0.36 + 0.05 8.09 0.28 2.59 0.01 
P011453500301 58011.200 0.28 + 0.12 8.00 0.57 2.02 0.01 
P011453500401 58012.260 0.29 + 0.06 9.96 0.43 2.76 0.00 
P011453500501 58013.255 0.5140.17 12.02 + 0.56 3.34 0.01 
P011453500601 58014.117 0.90 0.43 7.07 + 1.17 3.34 + 0.01 


P011453500901 58015.974 1.75 + 0.14 六 六 9.19 + 0.02 
P011453500902 58017.250 1.68 + 0.15 xk 9.36 0.01 
P011453500903 58017.389 1.85 + 0.35 ex 8.80 + 0.04 
GX 339-4 
P0304024015 59274.754 0.29 + 0.02 ex 1.89 + 0.03 
P0304024035 59295.758 0.04 + 0.03 ex 0.64 + 0.01 
MAXI J1348-630 
P0214002002 58511.314 0.08 + 0.02 6.39 + 0.26 0.4 + 0.01 
P0214002003 58512.308 0.11 + 0.01 7.07 + 0.24 0.53 + 0.01 


P0214002004: 01-08? 58513.303 0.15 0.01 8.95 0.25 0.76 0.01 
P0214002004: 09-19 58514.411 0.2 + 0.01 11.84 + 0.31 0.88 + 0.01 


P0214002005 58516.818 0.3 + 0.03 14.4 + 0.99 1.32 + 0.02 
P0214002006 58518.012 0.95 + 0.07 4.08 0.07 1.38 + 0.01 
P0214002007 58519.145 0.61 + 0.14 7.68 + 0.20 3.6 + 0.01 
P0214002008 58519.804 0.2 + 0.04 6.69 + 0.17 3.42 + 0.01 


MAXI J1631-479 
P0214003003 58527.1 x% 19.37 + 2.77 6.29 + 0.08 


P0214003008 58534.787 六 六 20.19+ 3.749 6.27 土 0.08 


GRS 1915 十 105 


P0101330006 58216.956 0.45 + 0.02 六 六 3.66 0.04 
P0101330007 58218.216 1.03 + 0.08 ex 3.96 + 0.55 
P0101330009 58220.272 0.32 + 0.03 ex 3.49 + 0.03 


a ** represent the PDS does not have this noise component in fitting. 


> The number behind the colon represents the exposure ID in this observation. 


25-6 


65 卷 


一 of PO11453500902 ME 
2 MJD:58017.25 i 
z 
3 
ae | 
2.5 | 
x 0.0 HT h, i i ual Ni 
K ataa NA 
-20 -15 -1.0 -0.5 00 05 10 I5 
v/Hz 
a 6 P011453500501 ME 
A MJD:58013.255 
E 
e | 
L o H 4 aptat ly 有 Ad 
tt pp | W | 
o is io òs 00 05 10 15 
v/Hz 
图 1 两 个 分 别 用 XSPEC 模 型 bknpow (上 ) 和 bkn2pow (下 ) 拟 合 噪 


声 的 PDS 


Fig.1 Two PDS in which the noise is fitted with the XSPEC 
model of bknpow (top) and bkn2pow (bottom), respectively. 


结果 
基于 上 一 章 的 数据 处 理 , 我 们 分 别 计 算 每 个 源 
在 3 个 探测 器 能 段 的 PDS, 拟 合 结果 列 于 表 2-4. 把 
MAXI J1535-571、 GX 339-4, MAXI J1348-630、 
MAXI J1631-479 和 GRS 1915+105 5 个 源 不 同 能 

的 截断 频率 和 QPO 频 率 , 也 即 久 与 yopo 的 关系 合 
并 于 图 2. 图 2 中 也 选取 了 Wijnands 等 回 的 结果 中 黑 
洞 的 数据 , 方便 与 我 们 的 结果 比较 . 从 图 2 可 以 看 到 
慧眼 -HXMT 的 3 个 能 段 和 Wijnands 的 数据 的 比较 
并 验证 W-K 关 系 , HXMT 观 测 到 的 数据 点 分 布 在 拟 
合 线 周围 , 离散 度 有 一 些 增加 . 特别 是 HE, 其 数据 
点 向 下 偏离 Wijnands 等 回 的 数据 . 从 图 2 中 可 以 看 


3 


马 宾 远 等 : 慧眼 -HXMT 观 测 到 的 


星系 统 中 的 W-K 关 系 


黑洞 双 


出 X 射 线 双星 的 与 yopo 之 间 存 在 窜 律 关系 , 对 于 
所 研究 的 观测 数据 , 我 们 利用 vopo = nvë (2 和 大 都 
是 此 窜 律 函数 的 参数 ) 拟 合 不 同 能 段 的 W-K 关 系 曲 
线 . 考虑 截断 频率 和 QPO 频 率 的 误差 , 拟 合 出 的 
参数 分 别 为 LE: n = 3.72 + 0.19, k = 0.69 + 0.07; 
ME: n = 5.61 + 0.50, k = 0.95 + 0.15; HE: 7 = 
6.08 + 0.60, k = 0.94 + 0.12; Wijnands 等 加 的 历史 
数据 拟 合 参 数 : n = 5.33 + 0.17, k = 0.88 + 0.07. 
拟 合 结果 也 显示 在 图 2. 定性 上 , 慧眼 -HXMT 观 测 
数据 和 原本 的 W-K 关 系 基 本 一 致 . 定量 上 , W-K 关 
系 离散 性 的 增加 以 及 拟 合 参数 的 差别 可 能 和 观 
测 数据 信 噪 比 有 关 , 例如 , 慧眼 -HXMT 3 个 探测 
器 较 短 的 曝光 时 间 ， 以 及 相 较 于 Wijnands 等 加 使 
用 的 探测 器 (例如 RXTE, Rossi X-ray Timing Ex- 
plorer) 更 小 的 低能 和 中 能 探测 器 的 有 效 面积 ; 另 一 
个 可 能 性 是 在 慧眼 -HXMT 观 测 数据 的 PDS 上 并 不 
能 很 好 地 区 分 不 同 成 分 , 例如 , MAXI J1820+070 
在 LHS 的 PDS 上 的 噪声 成 分 需要 分 成 好 几 部 分 来 
分 析 B031; 另外 , QPO 频 率 和 截断 频率 的 关系 可 
能 存在 本 征 的 离散 性 , 例如 , W-K 关 系 的 数据 点 
在 atoll 源 的 样本 中 离散 性 增加 , 而 Z 源 则 完全 偏 
离 W-K 关 系 [6, 3233], 因此 , 尚 需要 更 多 源 、 更 好 质 
量 的 观测 数据 验证 W-K 关 系 的 普 适 性 . 

为 了 演示 截断 频率 与 能 量 的 相关 性 , 图 3 展示 
了 MAXI J1535-571 源 宽带 噪声 的 高 频 截 断 频 率 
低频 截断 频率 加 随 能 量 的 变化 . 结果 表明 ， 
ww 对 能 量 并 没有 显著 的 依赖 性 , 这 和 在 MAXI 
J1820 十 070 所 发 现 现象 不 是 很 一 致 ,可 能 和 观测 
数据 信 噪 比 相关 . 与 之 相反 , ww 在 20 keV 以 上 随 
着 能 量 增加 有 增加 的 趋势 . 需要 说 明 的 是 ,慧眼 - 
HXMT 不 同 探 测 器 的 能 量 范 围 具有 重合 . 在 这 些 重 
BWAR, 不 同 探测 器 由 于 有 效 面积 和 本 底 水 平 的 
不 同 而 得 到 不 同 信 噪 比 的 探测 结果 , 我 们 在 此 处 以 
及 后 面 的 分 析 中 仅 保 留 信 噪 比 高 的 结果 . 

为 展示 快速 时 变 幅 的 变化 , 选用 典型 
源 MAXI J1535-571 的 录 声 成 分 的 rms 随 
能 量变 化 示 于 图 4. 在 低能 段 时 宽带 噪声 rms 保 持 接 
近 常 数 , 但 在 高 能 段 有 随 能 量 增加 而 下 降 的 趋势 ， 
此 结果 和 一 般 黑洞 LMXB 的 一 致 39. 
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图 2 从 LE、ME 和 HE 观察 到 的 5 个 BHB (红色 符号 


Fig.2 Correlation between low break frequency and QPO frequency of 5 BHBs (red symbols 


) 的 低频 截断 频率 与 QPO 频 率 的 关系 . 黑色 数据 点 取 


实 线 , 以 便 进行 比较 . 


Wijnands 等 [可 ,对 应 的 线性 拟 合 均 为 


from the observation of LE, ME 


and HE. The data (black points) are from Wijnands et al.[61. The corresponding linear fittings are all solid lines shown for 
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随 能 量 而 增加 Yang 等 B1 利 用 HXMT 
观测 数据 研究 了 MAXI J1820+070 的 PDS， 
高 能 段 的 PDS 的 截断 频率 也 是 光子 能 量 的 
频率 随 能 量 增加 . 我 们 对 5 个 黑洞 PDS 的 而 
明 , 虽然 离散 度 较 大 且 截 断 频率 在 每 个 探 
测 器 中 都 出 现 了 随 着 能 量 提升 而 增加 的 现象 , 但 在 
不 同 能 段 的 线性 关系 近似 保持 不 变 且 与 Wijnands 
等 四 一 致 . 这 表明 , 在 我 们 研究 的 观测 数据 信 噪 比 
水 平和 探测 器 精度 上 , 截断 频率 虽然 有 能 量 依赖 性 
并 没有 显著 改变 W-K 关 系 . 
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3 MAXI J1535-571 中 xb 和 上 与 能 量 的 关系 (观测 号 : 
P0114535003) 


Fig.3 The correlation between vp, Vn and energy in MAXI 
J1535-571 (ObsID: P0114535003) 
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4 MAXI J1535-571 中 以 下 的 低频 噪声 成 分 的 rms 谱 (观测 号 : 
P0114535003) 


Fig.4 The rms spectrum of low frequency noise component 
below vp in MAXI J1535-571 (ObsID: P0114535003) 


眼 -HXMT 观 测 到 的 黑洞 双 


星系 统 中 的 W-K 关 系 


不 同类 型 致密 天 体 XRB 符 合 相 同 的 W-K 关 系 ， 
表明 在 这 些 不 同类 型 XRB 中 , 吸 积 盘 的 成 分 相似 ， 
且 非 周期 性 的 快速 时 变 行为 可 能 具有 相同 的 物理 
机 制 . 对 宽带 噪声 和 QPO 变 化 幅度 的 研究 发 现 了 在 
低 硬 态 噪声 rms 高 达 50% 的 情况 , 表明 噪声 不 可 能 
起 源 于 标准 的 薄 吸 积 盘 B85, 噪声 的 延迟 性 质 也 部 
分 说 明了 这 一 点 B3. 这 些 快速 时 变 行为 可 能 产生 
于 非 热 辐射 成 分 的 内 吸 积 流 ( 冕 ). 对 于 噪声 的 形成 
机 制 , 人 们 建立 了 不 同 模型 861. 现在 的 主导 模型 则 
是 物质 吸 积 率 扰 动 传播 模型 , 在 此 模型 中 质量 吸 积 
率 波 动 由 粘 灌 变 化 所 传递 , 远离 黑洞 的 扰动 调制 趋 
于 稳定 并 接近 黑洞 的 吸 积 率 . 因此 , 内 吸 积 流 中 不 
同位 置 处 的 扰动 登 加 形成 我 们 观测 到 的 PDS. 在 这 
个 模型 中 , 两 个 截断 频率 分 别 对 应 不 同 的 吸 积 流 特 
FE 半径 , 低频 截断 频率 ww 和 吸 积 流 外 半径 rout (也 就 
是 吸 积 盘 内 半径 ) 相 关 , 而 高 频 截断 频率 思 则 对 应 
BAI ARAB Erin, 大 于 办 的 噪声 成 分 的 幅度 会 迅速 降 
低 ( 如 图 1). 

两 个 截断 可 能 符合 冕 上 的 同一 个 粘 滞 频 率 公 
式 B738]. 我 们 对 不 同 XRB 的 研究 表明 , 冕 的 分 布 是 
会 变化 的 , 随 着 能 谱 拟 合 参数 中 温度 的 升 高 , BY 
外 半径 减 小 , 这 和 等 离子 体 园 随 着 半径 增加 , 温度 
下 降 的 认识 一 致 . 我 们 的 结果 同时 也 表明 , 冕 内 半 
径 7i, 在 演化 中 保持 稳定 , 这 可 能 表明 rn 接近 最 内 稳 
定 轨道 半径 (Rirsco)ls9. 
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源 , 另外 也 可 能 是 现在 的 望远镜 观测 数据 不 足以 区 
分 不 同 的 理论 模型 , 导致 不 同 能 谱 模 型 的 简 并 , 所 
以 得 到 了 不 同 的 吸 积 盘 结 构 的 ]. 综合 利用 X 射 线 源 
的 时 变 和 能 谱 性 质 成 为 进一步 研究 吸 积 盘 结构 的 
途径 B0 43-44]. 为 此 , 我 们 选择 数据 较 丰 富 的 MAXI 
J1535-571 黑 洞 暂 现 源 , 通过 时 变 和 能 谱 分 析 研 究 
吸 积 盘 结 构 及 其 演化 . 

Kong 等 2 和 Huang 等 2 已 经 对 此 源 进 行 了 能 
谱 分 析 , 本 文 直接 借用 他 们 的 能 谱 参 数 计算 吸 积 盘 
内 半径 . 假设 吸 积 盘 与 热 内 流 稳定 衔接 , 那么 , 在 
吸 积 率 扰动 模型 中 , mm 应 与 热 内 流 外 半径 rou 也 就 
是 吸 积 盘 内 半径 反 相 关 . 图 5 显示 了 MAXI J1535- 
571 的 rou 与 加 的 相关 性 B3. 统计 上 , 图 5 中 截断 半 
径 和 截断 频率 的 确 呈 现 反 相关 关系 (其 中 不 同 颜色 
代表 了 不 同 探测 器 的 加 ), 支持 黑洞 暂 现 源 爆发 态 
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检验 
热 内 流 受 粘 滞 影 响 以 粘 滞 时 标 向 黑洞 运动 , 和 
粘 滞 时 标 联 系 的 粘 滞 频 率 由 下 式 表述 B7-39]: 


a0 


(1) 


其 中 mise 是 热 内 流 上 任意 半径 处 的 粘 清 频 率 ，a 为 
Hitt SB, 2 是 半 厚 度 (semithickness), h 为 热 内 流 
厚度 , 7 为 半径 , 为 粒子 的 开 普 勒 频率 . ,为 截断 
半径 也 就 是 ro 位置 上 的 粘 混 频率 . 由 于 (2)2 和 六 
都 是 半径 的 圭 律 形式 函数 , 仅仅 计算 同一 时 间 
内 吸 积 流 上 各 处 的 粘 滞 频 率 和 半径 相关 关系 时 ， 
(1) 式 可 以 简化 为 r = Bye, 其 中 B 和 6 是 这 个 窜 
律 关系 的 参数 . 对 吸 积 流 磁 致 不 稳定 性 的 模拟 显 


示 [85 , 热 内 流 的 w 和 “都 是 rout 的 函数 . 我 们 假设 rout 


Wise = Q ( 


Meee a IEA E OEN RRE i y BEE Ta. A BLE HA HEA IES ME 

增加 , 吸 积 盘 内 半径 缩小 的 物理 图 像 . 同时 我 们 也 是 近似 不 变 的 , 也 就 是 符合 相同 的 () 趟 ,选择 
能 看 到 , 不 同 的 办 可 以 对 应 很 相近 的 roue, 这 可 能 说 。 QPO 频 率 均 相近 且 时 间 上 连续 的 旧 光 号 数据 分 
明 能 谱 并 没有 很 好 反映 吸 积 盘 内 半径 的 演化 , 但 大 为 一 组 (1) 式 拟 合 图 5 中 的 数据 由 于 在 MAXI 
的 演化 趋势 还 是 很 清晰 的 . J1535-571 中 不 同 能 段 的 内 变化 相对 其 他 黑洞 来 说 
不 大 , 而 且 本 文中 的 截断 半径 由 能 谱 拟 合 (慧眼 - 
Z E HXMT 全 能 段 ) 得 到 , NSS ZS A E a e 
+ ME 比 , 就 把 所 有 能 段 的 数据 一 块 拟 合 , 我 们 把 3 个 探测 
103 + HE 器 的 数据 联合 拟 合 , 并 且 排 除了 截断 盘 可 能 已 经 到 
达 最 内 稳定 轨道 而 且 没 有 高 频 截断 的 最 右 侧 的 3 个 
£ 观测 号 . 这 样 我 们 可 以 将 图 5 中 的 数据 点 进行 拟 合 ， 
3 考虑 数据 点 在 X、Y 轴 上 的 误差 , 得 到 的 公式 参数 
10? AB = 289 + 88, 8 = 1.07 4 0.32. 借用 拟 合 公 式 
和 3 个 探测 器 的 高 频 截断 频率 ww, 可 以 计算 热 内 流 

FEAA 的 内 特征 半径 mw, 如 图 6 所 示 . 

e 50 如 结果 显示 ,内 吸 积 流 的 内 半径 变化 较 小 
vp/Hz (图 6). 为 此 , 我 们 取 图 6 中 3 个 探测 器 数据 计算 出 
图 5 MAXI J1535-571 中 低频 截断 频率 和 热 内 流 外 半径 rt 的 相 的 内 半径 的 平均 值 , 得 到 mm = 26.7 km. 按 前 人 
关 关 系 . 图 中 黑色 直线 是 (1) 式 对 数据 的 拟 合 ( 右 下 数据 未 参与 拟 合 , 原 的 论文 取 黑 洞 质量 为 10 个 太阳 质量 (Mo), Wrin ~ 
因 见 4.3 节 ). 0.89R,, Rs 是 黑洞 的 史 瓦 西 半 径 . 在 标准 吸 积 盘 


Fig.5 Correlation between the low break frequency (vp) and 

the outer radius (rout) of hot inner flow in MAXI J1535-571. 

The black line is the fit of data in Fig. 2 with Eq. (1) (Some 
data in low right is not included in this fit, the reason is 


.3). 


shown in Sec. 4 
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模型 中 , 内 吸 积 流 经 过 Risco 后 无 法 稳定 存在 并 进 
入 黑洞 , 普遍 认为 高 频 截断 频率 产生 的 位 置 也 接 
近 ISCO. 对 于 极端 Kerr 黑 洞 ， Risco ~ Rg. 本 文 获 
得 的 mn ~~ Risco lo 民 Rs 的 结果 和 内 吸 积 流 模 型 以 及 
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黑洞 自 旋 a = 0.993 一 致 , 表明 本 文 给 出 的 mm 具有 ar 
一 定 合理 性 . + ME 
102 + HE 
{ LE 
40 | | ME = 
1 HE S 
G 
35 > 
at] | | 
= 
= 25 10° 
02 0.4 0.6 08 1.0 
20 | Vp/Hz 
9 47 MAXI J1535-571 的 W-K 关 系 . 黑色 线 是 根据 Ingram 等 8 引 的 
58009 58010 58011 58012 58013 58014 理论 计算 的 W-K 关 系 . 
MJ 
Fig.7 The W-K correlation of MAXI J1535-571. The black 
图 6 和 (1) 式 计算 出 的 图 5 中 参与 拟 合 的 观测 号 的 ri line is the theoretical W-K correlation from Ingram et al.[35]. 


Fig.6 The inner radius rin calculated by v» and Eq. (1) in 
observations fitted in Fig. 5 


借 由 上 面 的 简化 公式 , 可 以 推论 出 内 吸 积 流 的 
内 半径 近乎 不 变 , 根据 计算 出 的 rm 和 其 他 的 黑洞 物 
里 参数 , 可 以 用 公式 分 别 计算 截断 频率 和 QPO 频 
率 . 在 L-T 相 对 论 进 动 模型 中 [39, QPO 频 率 由 下 式 
决定 : 


5720) all = (rin/fou)? t] c 
PO ’ 
me AHA) ESE — (rin rout) 4] Be 

(2) 
其 中 , a 为 黑洞 自 旋 因子 , 热 内 流 的 面 密度 分 布 符合 


TIE 
内 半 


因子 为 6 的 肾 律 函数 , rm 和 rout 分 别 为 热 内 流 的 


的 热 


径 和 外 半径 . 按照 Ingram 等 人 的 模型 855, 黑洞 


内 流 面 密度 不 随 着 


径 而 改变 , WE = 0, 假 


定 热 内 流 内 半径 min 


固定 了 


26.7 km, QPO 频 率 vopo 


就 可 以 当做 rout 的 函数 来 计算 , 截断 频率 也 同样 , 因 
此 我 们 可 以 画 出 理论 上 的 W-K 关 系 . BRAT AI, 
其 中 线 是 理论 值 构成 , 数据 点 是 前 面 图 中 给 出 的 
实际 观测 数据 , 右上 角 的 点 没有 参与 图 5 中 的 拟 合 ， 
因此 得 到 的 理论 曲线 偏离 了 它们 . 其 他 数据 点 与 
LE 论 值 的 差异 , 除了 观测 数据 分 析 结 果 的 粗糙 , 也 
和 公式 参数 过 于 简单 的 设置 有 关 [9. 但 趋势 性 的 
结果 表明 截断 盘 演化 模型 可 以 定性 解释 QPO 产 生 


和 W-K 关 系 . 
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5 结论 

本 论文 利用 慧眼 -HXMT 对 5 个 黑洞 候选 体 的 
观测 数据 , 在 2-120 keV 能 段 内 研究 了 W-K 关 系 及 
其 能 量 依赖 性 . 在 慧眼 -HXMT 3 个 探测 器 的 能 量 
BA, 功率 密度 谱 的 低频 截断 频率 没有 显著 的 能 
依赖 性 , 这 和 慧眼 -HXMT 以 及 NICER 对 MAXI 


FEL 


J1820+070 的 观测 结果 较为 一 致 . 慧眼 -HXMT 所 
给 出 的 W-K 关 系 与 Wijnands 等 g 在 3-20 keV 所 展 


示 的 结果 一 致 . 为 研究 吸 积 盘 结构 , 我 们 将 W-K 关 
系 应 用 到 MAXI J1535-571 暂 现 源 中 并 利用 能 谱 分 
析 的 成 果 , 结果 显示 MAXI J1535-571 系 统 与 截断 
吸 积 盘 的 模型 一 致 . 我 们 检验 了 QPO 与 吸 积 盘 结 构 
的 关系 , 发 现 QPO 频 率 与 热 内 流 的 相对 论 进 动 模型 


一 致 . 为 进一步 探索 X 射 线 双 星 吸 积 盘 的 结构 , 需 
要 更 多 的 观测 样本 . 
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The W-K Relation in Black Hole X-ray Binaries from 
Insight-HXMT Observation 


MA Bin-yuan'? YANG Zi-xu!? LIAO Jin-yuan? QU Jin-lu!? 
(1 School of Astronomy and Space Science, University of Chinese Academy of Sciences, Beijing 100049) 
(2 Institute of High Energy Physics, Chinese Academy of Sciences, Beijing 100049) 


Asstract Through decades of observations and studies, the characteristics of black hole X-ray binaries 
(XRBs) have been gradually revealed. However, the structure of the accretion disk remains uncertain. In 
the power density spectrum of the BHXRB (Black Hole X-ray Binaries), there is a correlation between 
the break frequency and the quasi-periodic oscillation (W-K relation), and it can constrain the current 
model of accretion disk. The W-K relation in black hole XRBs is studied based on the observations of five 
black hole XRBs by Insight-HXMT (Hard X-ray Modulation Telescope). The result shows that the W-K 
relation is valid in all the three telescopes with different energy bands. Moreover, a correlation between 
the break frequency and the inner radius of the accretion disk is found for MAXI J1535-571, which is 
self-consistent with the model of a truncated accretion disk. If the observed power density spectra result 
from the propagation of perturbation in mass accretion rate, it can be further inferred that the inner 
radius of the accretion disk is close to the innermost stable circular orbit and the black hole might be a 
high-spin system. 


Key words X-ray: binary, relativistic processes, accretion disk, radiation mechanisms: general 
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